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Abstract: We present an exact model of the detection statistics of 
a probabilistic source of photon pairs from which a fast, simple and 
precise method to measure the source's brightness and photon channel 
transmissions is demonstrated. We measure such properties for a source 
based on spontaneous parametric downconversion in a periodically poled 
LiNbC>3 crystal producing pairs at 810 and 1550 nm wavelengths. We 
further validate the model by comparing the predicted and measured values 
for the g( 2 \Q) of a heralded single photon source over a wide range of 
the brightness. Our model is of particular use for monitoring and tuning 
the brightness on demand as required for various quantum communication 
applications. We comment on its applicability to sources involving spectral 
and/or spatial filtering. 
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1. Introduction 

Sources of photon pairs are an essential building block in implementations of Quantum Com- 
munication protocols. Examples of such are Quantum Key Distribution (QKD), enabling un- 
conditional security in the exchange of confidential messages [1,2], or Quantum Repeaters, 
needed to break the distance barrier of QKD [3,4]. Photon pairs obtained from Spontaneous 
Parametric Downconversion (SPDC) [5] or Spontaneous Four- Wave-Mixing (SFWM) [6] in 
nonlinear materials, or from atomic ensembles [7, 8], can be used to generate either entangled 
photons by a careful arrangement of two downconversion paths [9], or to create a single pho- 
ton source by announcing the presence of one photon through the detection of the other one (a 
Heralded Single Photon Source, or HSPS) [10]. 

All aforementioned sources are of probabilistic nature, i.e. the number of emitted pairs per 
time unit follows a given statistical distribution such as a Poissonian or thermal distribution. In 
most applications, it is beneficial or even essential to know the mean number of photon pairs 
jj. emitted per unit of time, a quantity that we shall refer to as the brightness. For entanglement 
based QKD, Ma et al. have shown that both the key generation rate and the maximum distance 
over which a secret key can be established can be maximized by properly tuning the bright- 
ness [11]. Another example is the security of QKD based on HSPS, which relies on the ability 
of the sender to assess the photon statistics in a precise way [12-15]. Also, de Riedmatten et al. 
have shown that the visibility in Bell-state measurements, which is a key element of proposed 



quantum repeaters, crucially depends on the brightness [16]. 

Assessing the brightness of a source of photon pairs is a non-trivial task when limited to lossy 
channels and non photon-number-resolving detectors. This problem can be solved provided one 
knows the exact value of the total transmission of all photon channels. However, evaluating 
the loss associated with coupling a single photon from free-space to a singlemode fibre is not 
simple. One technique requires mode-matching a probe laser to the single photon mode, but this 
can be imprecise and unpractical (see [17-19] as examples). The brightness can also be inferred 
from measurements of the second-order autocorrelation function, g^ (0) [20]. However, as the 
time required for g^ (0) measurements depends on three-fold coincidence detection stemming 
from two simultaneously generated pairs, such measurements are time consuming to implement 
(see [8] and [21] as examples). Therefore, a method from which the brightness and the losses 
of the transmission lines can be determined with precision, speed and simplicity is necessary. 

In this work, we show how one can assess the brightness and the photon channel transmis- 
sions of a source of photon pairs by solely measuring single and two-fold coincidence detec- 
tions stemming from photons belonging to one pair. This makes this method very fast and 
efficient. In section|2l we present a model describing the detection statistics of any probabilistic 
source of photon pairs, and we show how the brightness and the losses of the channels can be 
assessed. Then, in section [3] we describe an implementation of the proposed method and use 
it to predict the value of the autocorrelation function, g( 2 '(0), of a HSPS for a wide range of 
brightness values. We confirm the model by the direct measurement of the predicted g^ 2 \0) 
values. Finally, in section|4] we discuss the limits of our model when the generated photons are 
spectrally correlated. 

2. An exact model of the detection statistics 

2.1. Description of the model 

To assess the properties of a source of photon pairs, we developed an exact model of the detec- 
tion statistics of the experimental setup detailed in Fig.[T] 



Fig. 1. The sources of photon pairs we consider comprise all probabilistic sources, includ- 
ing those based on nonlinear crystals, optical fibres or atomic ensembles. The distribution 
of the number of produced photon pairs per measurement time window can be given by 
any distribution such as Poissonian or thermal and is assumed to be known in advance. The 
pairs are deterministically separated, potentially by a dichroic beamsplitter in the case of 
collinear generation with non-degenerate wavelengths, or by non-collinear generation, into 
two separate channels. Each beam is filtered to remove all pump light and then the pairs 
are coupled into optical fibres. One beam is split again at a 50/50 beamsplitter before the 
photons are detected by non-photon number resolving single photon detectors and 



To model the detection statistics of this experimental setup we construct a column vector P, 
as shown in Eq. (Q]), which describes the joint state of the detectors: 
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P = ( Pabh Pabh Pabh Pabh Pabh Pabh Pabh Pabh ) . 



(1) 



Each element of P describes the probability that a set of detectors clicked or not per measure- 
ment time window, which is defined as the elementary observation time (e.g. a short time win- 
dow centered on one pump pulse; see later) for which detections are considered for statistical 
analysis. For example, p A M is the probability that detector Da clicked, during the measurement 
time window, and Dh and Dg did not. The goal is to determine how this vector, initially in state 
Po = ( 1 o ... o) T , is affected by single and multiple photon pair emissions as well as detector 
dark counts during one measurement time window. First, we describe the interaction of one 
photon pair with the detectors using the following transition matrix: 
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Each element of describes the probability for a pair to cause a transition of the three 
detectors. Each term is written as a function of T]h, t?a and tjb which are the overall trans- 
missions of each channel, from the photon pair source to Da and Dg respectively, in- 
cluding all optical losses, fibre coupling losses, detector inefficiencies, and the 50/50 beam- 
splitter. For example, M^(l,l) is the probability for the system to make a transition from 
ABH to ABH (i.e. to remain in the state where no detectors clicked), which must equal 

Pah + Pbh = (l-r}H-ilA + 1)aT)h) + (1 -r\H- i?b + J?bT?h) = 1 - t\h +_(rj A + J) B ){r\H - 1). 
Similarly, (2,1) is the probability to make a transition from ABH to ABH (i.e. no detectors 
clicked before and, provided one photon pair arrives, only Da clicks), which equals tj^(1 —"Hh)- 
All the upper diagonal elements are equal to as photons cannot make detectors "unclick". 
The rest of the matrix is constructed following the same physical reasoning. Furthermore, to 
conserve the total probability, each column of sums to 1 . The result of one photon pair 
interacting with the detectors is thus given by M^Pq- 

Second, the evolution of the system when ; photon pairs are created during the measurement 
time window is described by (Mn)'Po, as losses and the beamsplitter choice for individual pairs 
in multi-pair emission are independent. 

In addition to the absorption of a photon, thermal excitations can also cause detector clicks. 
These dark counts can be taken into account by constructing another matrix M ( / c . Thus, the 
evolution resulting from dark counts and ; photon pairs is described by M^ c (Mn) ! Po. Noting 
the dark count probabilities per measurement time window as du, cIa and ds, we get 
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Thus, when an unknown number of photon pairs are incident, it is possible to calculate the 
final vector P through 



p = Y,PiMdc(M n yp 



(4) 



where p\ is the probability to create i photon pairs per measurement time window. Provided 
that the probability distribution for p, is known this equation holds for all distributions, such as 
Poissonian, thermal or any distributions between the two [22]. Note that all matrices commute 



so the order in which they are applied does not matter. The construction of the matrices ensures 
that all elements of P are bounded individually between and 1 and that the elements of P sum 
to 1, i.e. the total probability is conserved. We note that the model is exact and there are no 
approximations. 

2.2. Determining the photon channel transmissions 

We now show how one can determine precisely the values of fi, r\n, "Ha and "Hb by measuring 
single and two-fold coincidence detection probabilities stemming from single pairs only. For 
these measurements, we require that the pump power (or equivalently the brightness of the pho- 
ton pair source) is low enough so that multi-pair events are negligible: /?,■ <C p i for i > 1 . Exper- 
imentally, this can be verified by looking at how correlated detections on D A are to detections 
on Dh- To measure this, we first define pn to be the heralding probability, i.e. the probability 
for Dh to click independent of the other detectors, pu = p A jj H + Pabh + Pabh + Pabh, and the 
similar expressions for pah and pa- We then define a parameter G = pah/{phPa) quantifying 
the strength of the correlation between detections at Da and Dh- The model described by Eq. 
predicts that, for Poisson, thermal and in between distributions, the value of G equals one at a 
very low brightness, when the coincidences are dominated by dark counts and detections are 
uncorrected, and equals one again at high heralding probabilities, when the coincidence detec- 
tions stem mostly from multi-pair emissions and correlations are smeared out. In between, the 
value of G can go well above 1 and this is an indication that multi-pair emissions are negligible. 
As we show here, this allows to experimentally obtain an upper bound for jj, when proceeding 
as follows. First, the dark count probability per measurement time window for each detector 
is measured. Second, the pump power is lowered and the transmissions are optimized until a 
value of G significantly higher than 1 is measured. Third, a plot of G versus fi is produced 
numerically assuming that the fibre coupling is perfect and that there are no additional optical 
losses, thereby setting the values of T]h and rj A equal to the specified detection efficiency of 
the detectors. Finally, an upper bound for fi is obtained from this plot by identifying the largest 
value of ji that produces a value of G equal to the measured value. They key point is that, for 
a given ji and dark count probabilities, G is decreased towards 1 when the transmissions are 
decreased. Thus, using this method, the true value of jj. must be smaller than the upper bound 
as the transmissions are overestimated. This, in return, allows one to obtain a lower bound for 
the ratio, r = pi/pt>u of the probability of single pair emissions, p\, over the probability of 
multi-pair emissions, p,>i = 1 — po — p\. As an illustration, using T]h = 60% and t\a = 25%, 
corresponding to the detection efficiencies of our detectors, and using their respective measured 
dark count probabilities (see section O, we produced the solid line shown on Fig. [2] where we 
assumed a Poisson distribution, p\ = exp(— ju)/i ! //!. 

Once the pump power is properly set and the lower bound on r is satisfactory, Eq. (0|i can be 
truncated to ; = 1 and one can show that the probability for Dh to click on a photon and not a 
dark count is given by p^ = (pn — dn)/ (1 — da). Similarly, we get p A = (pa — d A )j (1 — d A ) 
and the equivalent for p^ \ In the same way, we can get expressions for the coincidence prob- 
abilities pah and pbh- Then, using these expressions and an experimental data collection run 
with a heralding probability that guarantees negligible multi-pair events, one can solve for the 
four unknowns ji, rj#, T]a and V]Bi since the dark count probabilities can be measured directly. 
These unknowns can be calculated through equations ((5]) through (0. The equivalent set for Dg 
is constructed by replacing v\a and p A by v\b and pg \ respectively: 

Pah- PiPdA(l -d H ) -p A 1) d H {l-d A )-d A d H ... 
p [ A ] (\-d A )(\-dH) 
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Fig. 2. Correlation strength G versus brightness jl. The solid line corresponds to T]h = 60% 
and rj A =25%. It reaches a maximum value at very low fl and then sharply decreases 
to 1 for fi = (not visible). The meanings of the dotted and dashed lines are discussed in 
section[3]and[4]respectively. 
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Note that these predictions apply to any statistical distribution for which the multi-pair events 
can be neglected (for example, through the method described above). However, to determine the 
value of the brightness, one must have prior knowledge of the distribution and how to relate it 
to the measured value of p\. In the case of a Poissonian source, we have p\ = jxexp(—n) 
which can be solved numerically for jj.. The case of a thermal distribution is similar with 
p\ = (tanhju/coshju) 2 . 

Once the transmissions are precisely known, one can use Eq. (01 to find the brightness that 
corresponds to any measured heralding probability. This will then allow to predict the complete 
detection statistics vector P. 

2.3. Application to a HSPS 

The transmissions and the brightness, along with the knowledge of the pair distribution type, 
can be used to predict the second-order autocorrelation function of the heralded mode of a 
HSPS, g^ (0), for any desired heralding probability pn, in a Hanbury Brown and Twiss (HBT) 
experiment [23]. The distribution of the number of photons in that mode follows the distribution 
of the number of photon pairs created by the source except for a reduced vacuum component, 
po, due to the heralding. A g( 2 > (0) < 1, which is achievable with a HSPS, implies a nonclassical 
source (for a perfect single photon source g^ 2 \0) = 0). Alternatively, a g^ 2 '(0) > 1 describes a 
classical source (for Poissonian g^ 2 \0) = 1 and for thermal g^(0) = 2). To verify our model, 
we compare its predictions with a real measurement of the g( 2 \0). In this experiment, which 
can be seen as measuring a subset of Eq. dU, detectors D A and Dg are activated only when Dh 
clicks. The g( 2 '(0) is defined as 



(8) 



where Pab\h is tne probability that both and Dg click provided that Dh clicked, etc. 

For a specific heralding probability pn, we can directly measure g' 2 '(0) using the setup of 
Fig. [T]by keeping only the events where Dh clicked. On the other hand, the g^ 2 \0) can also 
be predicted for the same heralding probability using Eq. (@). The experimental results of this 
verification are presented in the next section. 

One interesting result can be derived from our model. Considering a Poissonian distribution 
at low brightness and assuming that dark counts are negligible, one can derive from Eq. 
that g( 2 \0) — /J. (2 - 7]h) [24]. This indicates that for a HSPS, the transmission to the heralding 
detector is a crucial parameter to optimize. 

It is important to note here that the 50/50 beamsplitter used in the setup of Fig. [TJ is not 
required to determine the brightness and the transmissions. Indeed, to assess the brightness and 
the transmissions only the detectors Dh and Da are necessary. In this work, the beamsplitter 
and Dg were added only to provide a way to verify the validity of the predictions through the 
g( 2 ' (0) measurement. Modifying the vector P and the matrices to accommodate for a setup with 
no beamsplitter and Dg is straightforward. 

3. Experimental results 

The experimental setup is shown in Fig. [3] A clocking signal triggers a pulsed laser diode (Pi- 



coQuant PicoTA) creating 50 ps pulses at 532 nm filtered to remove excess 1064 nm light. 
The pulses are focused onto a 1 cm long periodically poled LiNb03 crystal (PPLN) from 
Stratophase with a grating period of 7.05 |jm heated to 175.7 °C. Then, collinear spontaneous 
parametric downconversion to one or several photon pairs can occur, each pair consisting of one 
810 nm and one 1550 nm photon. A dichroic mirror DM separates the two wavelengths and, af- 
ter removing the 532 nm light with long-pass color filters, the photons are coupled into SMF28 
optical fibres. The 810 nm photons are sent towards a free-running Si single photon counting 
module D H (SPCM-AQR-14-FC, Perkin-Elmer) and the 1550 nm photons are detected by two 
gated InGaAs single photon detectors D^ and Dg (id201, IdQuantique) positioned right after a 
50/50 fibre beamsplitter. The detection statistics are recorded using a Time-Digital-Converter 
(TDC-GPX, ACAM), providing the time elapsed between a start pulse, given by the clock, and 
each stop pulse, corresponding to detections. The data is analyzed in real-time and the statistics 
are updated continuously until the end of each run. The coherence length / ( of the downcon- 
verted 810 nm photons was measured to be 180 jxm using a Michelson interferometer. From 
these measurements we calculated the bandwidth to be AAsio ~ 4nm and, based on energy 
conservation of the SPDC process, AA1550 ~ 15 nm. As the downconverted photons' coherence 
time, which equals l c /c = 0.54 ps, is much smaller than the pump pulse duration, which is 50 ps, 
we can confidently assume that our source of photon pairs follows Poissonian statistics [22]. 
The InGaAs detectors, which require gating to limit excess dark counts, were activated from 
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Fig. 3. Experimental setup. 



each clocking signal for a 5 ns measurement time window. Detections on the Si detector were 
considered valid only if they arrived within a 5 ns window centered on the clocking signal, as 
measured by the TDC. To determine the transmissions, the clocking signal triggering the laser 
and InGaAs detectors was set to 30 kHz. This low repetition rate ensured that saturation effects 
in the detection electronics and biases in the detection statistics from the InGaAs detectors' 
10 [is dead-time were avoided. We first measured the dark count probabilities to be = 2.87 x 
10~ 4 , dg = 3.84 x 10~ 4 and du = 2.5 x 10~ 7 per 5 ns. Next, we lowered the pump power using 
neutral density filters in order to increase the correlation strength between Dh and Da to a value 
of G = 20.6 ± 1.0, corresponding to a heralding probability of 0.287 ± 0.001%. Intersecting 
this value with the solid line of Fig. [2] gives a upper bound of jU < 0.0480 ± 0.0013, yielding 
r > 41 .0 ± 2.2, which we considered sufficiently high to continue. Next we measured the single 
and coincidence detection probabilities from which we obtained the following values: T]h 
0.1212 ■ 0.0031. n A =0.0145 ±0.0005, J] B = 0.0162 ±0.0005 and /i = 0.02375 ±0.00016, 
corresponding to r = 83.5 ±0.6. The G curve corresponding to these values is plotted as the 
dotted line on Fig. [2] and the predicted value of G at jj. = 0.02375 is 23.9 ±0.5, which is close 
to the measured value of 20.6. 

Using these values together with Eq. ©, we produced a plot of the predicted Pab\h^ Pa\h an d 
Pb\h f° r a wide range of the brightness (and consequently, of the heralding probability). These 
predictions are compared to the measured values on Fig. Ufa) and@tb). Next we compared the 
predicted and measured (0) as shown on Fig. [2a). On the same figure we plotted the value 
of the brightness fi corresponding to each heralding probability. In all cases, the agreement 
between the predicted and measured values is excellent. We note that for these measurements, 
the repetition rate was increased to 5 MHz and the InGaAs detectors were activated for 5 ns 
only when the Si detector clicked synchronously (within a 5 ns window) with the pump, as 
required for g^ 2 '(0) measurements in the HBT setup. This resulted in an average detection 
rate of 30 kHz, with randomly distributed time differences, for the InGaAs detectors and was 
thus sufficient to ensure that saturation effects in the detection electronics was not an issue. 
However, to eliminate the effect of the dead-time on the detection statistics, we considered only 
the events where both InGaAs detectors were ready to detect photons, as provided by the "gate 
out" electrical signals of these detectors. 
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Fig. 4. (a) Predicted (solid lines) and measured (points) conditional detection probabilities 
Pa\h atM1 Pb\h- (b) Predicted (solid line) and measured (points) conditional coincidence 
probability Pab\h ■ The dashed lines on both plots are the one standard deviation uncertainty 
bounds on the predicted values which were generated using the uncertainty bounds on the 
measured transmissions. 



Our method drastically reduces the time needed to characterize the source as measurements 
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Fig. 5. (a) Predicted autocorrelation g' 2 ' (0) for Poissonian (solid line) and thermal (dotted 
line) distributions, measured values (points), and the corresponding brightness fi (dash- 
dotted line). The measured data agrees very well with the Poissonian distribution. The 
dashed lines are the one standard deviation uncertainty bounds on the predicted values, 
(b) As the heralding probability reaches the noise level of Dh (dashed line), the model 
correctly predicts that the g^ 2 '(0) approaches one, as uncorrected dark counts begin to 
dominate over photon clicks. 

of single and two-fold coincidence detections at a low heralding probability are sufficient to 
determine the transmissions. These can then be used to predict the brightness and the g^ (0) of 
a HSPS for any heralding probability. In contrast, a single direct measurement of the g^ (0) at 
a given heralding probability requires three-fold coincidence detections stemming from multi- 
pair emissions, which are less likely to happen. In our experiment, at a heralding probability 
of 0.287%, two-fold coincidences were approximately 700 times more likely than three-fold 
coincidences. Consequently, a direct g^ 2 \0) measurement required much more time. 

4. Effect of spectral and spatial correlations 

The model proposed in section[2]may not apply directly to a source when spectral and/or spatial 
correlations exist between the photons of each pair and when the transmission of each chan- 
nel are frequency and/or spatially selective. For instance, Bell state measurements generally 
requires filtering of photons that are spectrally correlated [25]. Let's suppose the spectral fil- 
tering applied on each spatially separated photons is performed with two separate filters that 
both need to be aligned on the photons spectra. Due to energy correlation, transmission of one 
photon through a spectral filter determines the spectrum of the other photon [26]. If the filter- 
ing of the second photon does not match its now modified spectrum, the coincidence detection 
probability is reduced. In this case, given that one photon pair was created, we can still write 
the probability to get a detection at Dh to be 7]h and at Da to be 7]a- However, the probability 
to get a coincidence is lowered to cv\h1]a, where < c < 1. The upper bound is reached when 
either the photons' spectra are uncorrected before the filtering (see [27]), or when the selected 
spectra satisfy the energy conservation conditions perfectly. 



When this situation applies, the detection matrix can be re-written as follows: 
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The presence of spectral correlations affects the predictions of the proposed model but as we 
show here, the consequences are minimal. First of all, for a given brightness, a value of c < 1 
lowers the measured value of G towards 1 . Therefore, when assessing if multi-pair events are 
negligible, one can assume that c = 1 and the upper bound on /i, along with the lower bound 
on r, are still valid. As an illustration, reducing c from 1, corresponding to the solid line of 
Fig. |2] down to 0.5, corresponding to the dashed line on the same figure, lowers the curve to- 
wards G = 1 . Also, performing the above analysis to calculate the transmissions, while neglect- 
ing dark counts, we can show that the obtained solutions are rj' H = ct\h, r\' A = cT]a, r\' B = crjg 
and fi' — fi /c, and that the value of c cannot be assessed directly. Therefore, if one is unaware 
of c, the analysis performed underestimates the transmissions by a factor of c and overestimates 
the brightness by a factor of 1/c. However, we can show by direct calculation that the predicted 
probability vector P, and consequently the predicted g^ (0) for a given heralding rate, are inde- 
pendent of c. Therefore, in our results obtained in section[3] we may have slightly overestimated 
the values of the brightness but we still predicted the correct value for the measured g^ (0). 

For QKD with a HSPS, one can assess the security of the protocol against PNS attacks by 
knowing the value of ji. The analysis we propose gives an overestimated value for jj. when 
there are reasons to expect that spectral correlations and unmatched filtering are present. This 
is not detrimental to the security of QKD, as a sender unaware of this will, in the worst case, 
only overestimate the information available to an eavesdropper and shorten the key more than 
necessary through privacy amplification [28,29]. 

5. Conclusion 

We developed a model exactly describing the detection statistics of a probabilistic source of 
photon pairs. From this model, we outlined a method from which the transmission of each 
photon channel, as well as the source's brightness, can be determined by measuring single 
and two-fold coincidence detection probabilities stemming from photons belonging to one pair. 
Then, we experimentally confirmed the method by demonstrating that the measured g^ (0) of a 
HSPS can be correctly predicted for any heralding probability. This allows one to quickly tune 
the brightness on demand as required to optimize the performance of entangled QKD, to assess 
the security of HSPS-based QKD or to optimize quantum repeater error rates and distances, 
all in the context of fluctuating experimental conditions such as photon channel transmissions. 
Finally, we showed that our model correctly reproduces the detection statistics even if the pho- 
tons are spectrally and/or spatially correlated, and that this only leads to an overestimation of 
the brightness of the source. The simplicity of the proposed method makes it very attractive for 
the field of quantum communication in general. 
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